function P = RungeKuttaMethod( Fxy, y0, a, b, n)
f = inline(Fxy, 'x', 'y');
h = (b - a) / n;
P = struct('x', num2cell(a : h : b), 'y', num2cell(zeros(1,n + 1)));
P(1).y = y0;
for i = 2 : n + 1
    fk = f(P(i-1).x, P(i-1).y);
    P(i).y = P(i-1).y + h * (fk + f(P(i).x, P(i-1).y + h * fk)) / 2;
end;
end

